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Abstract 

In this note, we present a natural proof of a recent and surprising result 
of Gregory Berkolaiko interpreting the Courant nodal defect as a Morse 
index. This proof is inspired by a nice paper of Miroslav Fiedler published 
in 1975. 

1 Introduction 

The "nodal defect" of an eigenfunction of a Schrodinger operator is closely related 
to the difference between the upper bound on the number of nodal domains 
given by Courant 's Theorem and the number of nodal domains. In the recent 
paper [2], Gregory Berkolaiko proves a nice formula for the nodal defect of an 
eigenfunction of a Schrodinger operator on a finite graph in terms of the Morse 
index of the corresponding eigenvalue as a function of a magnetic deformation of 
the operator. His proof remains mysterious and rather indirect. In order to get 
a better understanding in view of possible generalizations, it is desirable to have 
a more direct approach. This is what we do here. 

After reviewing our notation, we summarize the main result and give an in- 
formal description of the proof in Section |3j The proof itself is implemented in 
Sections H] and [5] with an alternative view provided in Appendix |A] The contin- 
uous Schrodinger operator on a circle is considered in Appendix [B] and various 
special cases and further ideas are explored in other Appendices. 
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2 Notation 



Let G = (X, E) be a finite connected graph where X is the set of vertices and E 
the set of unoriented edges. We denote by {x, y} the edge linking the vertices x 
and y. We denote by E the set of oriented edges and by [x, y] the edge from x 
to y; the set E is a 2-fold cover of E. A 1-form a on G is a map E — y R such 
that = — a([x, ?/]) for all {x, y} G E. We denote by the vector 

space of dimension j^E of 1-forms on G. The operator d : 1R X — > is defined 

by df([x, y]) = f(y) — f(x). If Q is a non-degenerate, not necessarily positive, 
quadratic form on Q l (G), we denote by d* the adjoint of d where M. x carries the 
canonical Euclidean structure and Q}(G) is equipped with the symmetric inner 
product Q associated to Q. We have dimkerd* = j3 where j3 = 1 + j^E — j^X is 
the dimension of the space of cycles of G. We will show later that, in our context, 
we have the Hodge decomposition Q l (G) = dM. x © ker d* where both spaces are 
Q-orthogonal. 

Following jl], we denote by Oq the set of X x X real symmetric matrices 
H which satisfy h XjV < if {x,y} G E and h x>y = if {x,y} ^ E and x ^ y. 
Note that the diagonal entries of H are arbitrary. An element H of Oq is called 
a Schrddinger operator on the graph G. It will be useful to write the quadratic 
form associated to H as 

?i(/) = - E w/(^)-/(2/)) 2 + E^( x ) 2 ' 

{x,j/}G-B xSX 

with = h x , x + magnetic field on G is a map B : E U(l) 

defined by B([a;, y]) = e iax ' v where [x,y] h-> a X)!/ is a 1-form on G. We denote by 
Bg = e in ^ the manifold of magnetic fields on G. The magnetic Schrddinger 
operator Hb associated to H G Og and 5 = e* Q is defined by the quadratic form 

Mf) = -\ E ^ /(•'•) '/(//) 2 • E r --/v 

[x,2/]eE ^gx 
associated to a Hermitian form on C x . More explicitly, if / G C x , 

///(•'•) = WW + E h x , y e ia *«f(y) . (1) 

We fix i7 and we denote by 

Ai(fl) < A 2 (S) < < A n (S) < • • < A # x(fi) 

the eigenvalues of Hb- It will be important to notice that X n (B) = X n (B). 
Moreover, we have a gauge invariance: the operators Hb and Hb' with a 1 = a+df 
for some / G M x are unitarily equivalent. Hence they have the same eigenvalues. 
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This implies that, if Q 1 (G) = d'R x © kerd* (this is not always the case because 
Q is not positive), it is enough to consider 1— forms in the subspace kerd* of 
Q 1 (G) when studying the map A n : B — > X n (B). This holds in particular for 
investigations concerning the Hessian and the Morse index. 

3 Statement of Berkolaiko's magnetic Theorem 

Before stating the main result, we recall the 

Definition 1 The Morse index j(q) G NU {+00} of a quadratic form q on a real 
vector space E is defined by j(q) = sup F dimF where F is a subspace of E so 
that q\F\o is < 0. The nullity of q is the dimension of the kernel of q. 

The Morse index of a smooth real-valued function f defined on a smooth 
manifold M at a critical point Xq £ M (i.e. a point satisfying df(xo) = 0) is the 
Morse index of the Hessian of f , which is a canonically defined quadratic form 
on the tangent space T X0 M. The critical point Xq is called non-degenerate if the 
previous Hessian is non-degenerate. The nullity of the critical point x$ of f is 
the nullity of the Hessian of f at the point xq . 

The aim of this note is to prove the following nice results due to Berkolaiko 
[HE]: 

Theorem 1 Let G = (X, E) be a finite connected graph and (3 the dimension 
of the space of cycles of G. We suppose that the n-th eigenvalue \ n of H G Oq 
is simple. We assume moreover that an associated non-zero eigenfunction <p n 
satisfies <fr n (x) 7^ for all x G X . Then, the number v of edges along which <j) n 
changes sign satisfies n — 1 < u < n — 1 + /3. 

Moreover A n : B — > \ n (B) is smooth at B = 1 which is a critical point of A n 
and the nodal defect, 8 n — v — (n — 1) is the Morse index of A n at that point. 
If M is the manifold of dimension (3 of magnetic fields on G modulo the gauge 
transforms, the function [B] — > A n (B) has [B = 1] as a non-degenerate critical 
point. 

Remark 1 The previous results can be extended by replacing the critical point 
B = 1 by B x y = ±1 for all edges {x, y} G E. The number v is then the number 
of edges {x, y} G E satisfying B X: y(f) n (x)(f) n (y) < where <p n is the corresponding 
eigenfunction. 

Remark 2 The assumptions on H are satisfied for H in an open dense subset 
ofO G . 

The upper bound of v in the first part of Theorem [T] is related to Courant nodal 
Theorem (see [B] Section VI. 6) as follows: a nodal domain on a graph for the 
eigenfunction <p n is a connected component of the sub-graph G' of G obtained 
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by removing the edges along which <p n changes sign. Denoting by [i the number 
of nodal domains of n , the Courant Theorem for graphs (see [I], Theorem 2.4) 
asserts that \i < n; using Euler formula for the graph G' and because /i = bo(G'), 
the number of connected components of the graph G', we get also a lower bound 
(see P): 

Corollary 1 Under the assumptions of Theorem^ we have n — (3 < fi < n . 

Important warning: Without loss of generality, we can and WILL assume 
in the rest of this note that X n = A n (l) = 0. This implies that the Morse index 
of q± is n — 1 . 

In the course of the proof we will use a special choice of gauge in which we 
can compute the Hessian explicitly. More precisely, according to the classical 
perturbation formulae, 

^ = (0> H(j)) + 2(H(f), 0), 

where we assumed that A is at a critical point: A = 0. The first term is easy to 
calculate explicitly; for perturbation in the direction of the 1-form u it is 

Q{u) = o X/ a x,y u (l x ' y\) 2 with a *,y = -h x ,y4>n(x)4>n(y) = a y ,x ■ (2) 

E 

Considered as a quadratic form in oj, Q is already in the diagonal form. Its index 
is clearly the number of negative values among {— h Xty 4> n (x)(j) n (y)}, or, in other 
words, the number v of edges where n changes sign! 

We will present an explicit choice of gauge in which the second term vanishes. 
The condition for this is Hcf) = which, after explicit calculation, can be inter- 
preted as u G kerrf*, where d* is the conjugate of d with respect to the inner 
product induced by (J2J). 

Finally, we observe that the index of Q(u) has been computed to be v in the 
whole of Q 1 (G), whereas we should be restricting ourselves to our chosen gauge, 
u G ker d*. We will show that this restriction reduces the index precisely by n — 1. 
Indeed, the splitting Q 1 (G) = dWL x ©ker d* is orthogonal with respect to the form 
Q, therefore 

ind(Q) = ind (Q\<mx) + ind (Q\ keid *) ■ 

We establish that ind(Q|diRx) = n — 1 by relating the form Q on dM. x to the 
quadratic form qi around the point <p n . 

4 The quadratic form Q 

Lemma 1 The set of forms f — > (f(x) — f{y)) 2 where {x, y} G V2(X), the set 
of subsets with two elements of X , and f — > f(x) 2 with x G X is a basis of the 
set of quadratic forms on M. x . 
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Definition 2 A quadratic form q on M x is said of Laplace type if V/ G q(l, f) 
where q is the symmetric bi-linear form associated to q. 

Lemma 2 TTie set of forms f —> (f(x) — f(y)) 2 , {x,y} G V%{?C) is a basis of 
the space of quadratic forms of Laplace type. 

The form q± : f — > qi(4> n f), where 4> n f is the point- wise product of (f> n and /, 
is of Laplace type because 

?i(l,sO = {H<pn\<p n g) = (O|0„p) . 

Hence gi(l, g) = 0. 

Moreover, qi(f) = Q(df). Indeed, because of Lemma [2} it is enough to 
compare the coefficients of the basis forms / — > (f{x) — f(y)) 2 . The form / — > 
Q{df) is already expanded in this basis. To find the coefficient for the form 
/ ~~ Qi{f)i we observe that (because we know it is of Laplace type) the coefficient 
in question is minus the coefficient in front of the term f(x)f(y), divided by two. 
This evaluates to a XjV (see equation §2§). 

In fact, we will need to use Q(df,dg) = (H((j) n f)\(j) n g). 

Lemma 3 The Morse index of Q\dRX is equal to n — 1. 

It is a general fact that the Morse index of the quadratic form / — > Q(Af) is 
the same as the Morse index of the restriction of Q to the image of A. Hence, 
the Morse index of Q\dM.x is the Morse index of q~\ on M. x . Because / — > <p n f is a 
linear isomorphism, this index is equal to the index of q\ by Sylvester Theorem. 
Since A n = 0, the index of q\ is n — 1 by elementary spectral theory. 

Lemma 4 Let us denote by d* the adjoint of d where M. x is equipped with the 
canonical Euclidean structure and Q 1 (G) with the inner product associated to Q. 
The space fi 1 (G) splits as 

Q\G) = dR x ©kerd* 
(Hodge type splitting), and this decomposition is Q-orthogonal. 
More explicitly d* is given by 

d*u{x) = ^2a Xjy u)([y,x)) . 

If u — df satisfies d*u = 0, we have d*df = 0. Hence Q(df, dg) = for all 
g and {H((j) n f)\(j) n g)) = 0. Because A n is of multiplicity 1, this implies that / is 
constant and hence df = 0. So cflR x n ker d* = {0} and the conclusions follow. 

At this point, we know that the nodal defect is the Morse index of the re- 
striction of Q to the space ker<i* of dimension 0. The first part of the Theorem 
follows. 
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5 The magnetic Hessian 



We need one more fact to complete the proof: to identify the Hessian of A n on 
e ikerd b = 1 with the restriction of Q to kerd*. 

Let us denote by S C C x the set of unit vectors / normalized so that f(x ) 
is real and f(xo) > where Xq is chosen in X. 

Lemma 5 The point B = 1 is a critical point of A n . If <f) n (B) G S is the 
eigenf unction of Hb corresponding to the eigenvalue X n (B), the differential of 
B — > 4> n (B) vanishes at B = 1 on kerd*. 

The first property comes from the fact that A n (B) = A n (B). We can compute, 
for any variation e lta , t close to 0, of B = 1, Hb4>u + H<fi n — 0. The condition 
d*at = can be written as Yl y ~x ^ l x.y4'n{y) < ^x,y — for all x G X. From Equation 
([I]), this is equivalent to HB<t>n = 0. Hence H((j) n ) = and n = c0 n since 
X n is simple. From the normalization ||0 n (5)|| = 1, we get c G iM. and, since 
(f> n ( x o) ^ the number c is real. We deduce that <p n = 0. 

Lemma 6 The function F : S x e ikcrd * -> R rfe/ined 5y F(f,e ia ) = (H eia f\f) 
admits (<f) n , 0) as a critical point and the Hessian of (A n )| e ik Cr d* at the point B = 1 
is £/ie /orm Q. 

The differential of F with respect to / vanishes because / is an eigenfunction of H . 
The differential with respect to kerd* vanishes, because F(f,e ia ) = F(f,e~ ia ). 
The Hessian of F at (<f) n , 0) is well defined. Because the differential at B — 1 
of B — > <p n (B) vanishes on e lkcrd \ the Hessians of : B — > F(<f> n (B),B) and 
M n : B — > F((f> n (l),B) agree. A simple calculation of the Hessian of M n gives 
the result: 

[x,y]£E x ^ x 

- ^2 hx >v (^™( x ) 2 + ^{yf - 2 cos a X)V (j) n (x)<j) n (y)) + ^ V x \(p n (x)\ 2 . 

[x,y]£E x£X 

Computing the second derivative with respect to a at a = gives Hessian(M n ) = 
Q(a). 



A A pedestrian approach to the calculus of the 



Hessian of A n in Section [5 



We will derive a direct approach to the calculus of the second derivative of an 
eigenvalue which could be used directly in the proof of Lemma [61 Let t — > A(t) 
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be a C 2 curve defined near t = in the space of Hermitian matrices on a finite 
dimensional Hilbert space ("H, (.|.)). Let us assume that A(0) is an eigenvalue of 
A(0) of multiplicity one with a normalized eigenvector 0(0). Then, for t close to 
0, A(t) has a simple eigenvalue X(t) of multiplicity one which is a C 2 function of 
t. We can choose an associated eigenfunction 0(t) which is C 2 with respect to t. 
The following assertions give the values of the first and second derivatives of X(t) 
at t = 0: 

Proposition 1 Under the previous assumptions, we have 

A'(0) = <A'(O)0(O)|0(O)> , 

J/A'(0) = 0, we /iawe 

A"(0) = (A"(O)0(O)|0(O)} + 2<0'(O)|A'(O)0(O)> , 

w/iere 0'(O) is any solution of (A(0) - A(O))0'(O) = -A'(O)0(O). 
In particular, if A' (0) 0(0) = 0, 

A"(O) = (A"(O)0(O)|0(O)> . 

Proof. - 

We start with (A(£) — X(t))<f)(t) = where 0(t) is an eigenfunction 
of A(t) which depends in a C 2 way of t. Taking the first derivative, 
we get 

(A'(t) - A'(t))0(t) + (A(t) - \{t))4f(t) = . (3) 

Putting t = and taking the scalar product with 0(0), we get the 
formula for A'(0). Similarly, the t-derivative of Equation ([3]) is 

(A"(t)-A"(t))0(t)+2(A'(t)-A'(t))0'(t) + (A(t)-A(t))0"(t) = . (4) 

Pouting t — 0, taking the scalar product with 0(0) and using A'(0) = 
0, we get the result. 

□ 

We can apply this to A(t) := H e it a with a G ker<i* in order to get the Hessian 
of A n in Section [51 The condition A'(O)0(O) = is exactly d*a = 0! 

B Hill's operators 

In this Appendix, we will describe the case of a Schrodinger operator on the 
circle, also called the Hill's operator. This is the simplest continuous case, but it 
may be useful to do it with some details in order to try to extend the method to 
higher dimensional manifolds. 
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Eigenvalues and discriminant 

The Hill's operator is 

d 2 

H = -— + q(x) 
dx z 

where q : R — > R is a smooth, 1-periodic, function. The spectral theory of Hill's 
operators has been well studied; in particular, the inverse spectral theory for 
this operator allows to solve non-linear evolution equations, like the Korteweg-de 
Vries one. A presentation of the properties of Hill's operators is given in [TP] . 
The following facts are known: 

Theorem 2 If we denote by A • , j — 1, • • ■ the spectra of H acting on periodic 
(resp anti-periodic) functions of period 1, we have the inequalities 

Xt<K< A 2 - < A+ < A 3 + < • • • 

and the spectrum of H on L 2 (R) is then union of intervals, called the bands, 

[A+,A r ]U[A 2 -,A+]U[A3 + ,A 3 -]U--- . 

These statements are linked to the properties of the discriminant A (A): if yi(x, A) 
and 1/2 (x, A) are the normalized solutions of (H — \)y = whose Cauchy data 
are ?/i(0, A) = 1, y[(0,\) = 0, 1/2(0, A) = 0, y f 2 (0,\) = 1, the discriminant A is 
the entire function given by A(A) := yi(l, A) + y' 2 (l, A). The spectrum of H on 
L 2 (R) is the set of real A's so that |A(A)| < 2. The periodic (resp. anti-periodic) 
spectra are given by A(A) = 2 (resp. A(A) = —2). The function A(A) — 2 is a 
regular izat ion of n^Li(-^ — Ki) lY1 the spirit of [5]. It is proved in [10], Section II, 
that, if A+ is simple, A'(A+) 7^ and the sign of this derivative is that of (—1)™. 

Magnetic fields 

We will assume that A+ is equal to and is a simple eigenvalue of H acting on 
1-periodic functions. Up to gauge transform, every magnetic potential on the 
circle is a constant a. The bands are linked to the addition of a magnetic field 
as follows: the n-th band is the image of the circle U = {e ia \a G R} by the map 
A n where A n (e ia ) is the n-th eigenvalue of H a which is H acting on functions / 
so that f(x + 1) = e ia f(x). In particular, if n is even, A^ is a maximum of A n 
while if n is odd, A^ is a minimum of A n . This fits with Berkolaiko's formula 
because the (even!) number of zeros of the corresponding periodic eigenfunction 
4> n is n = (n — 1) + 1 if n is even and n — 1 = {n — 1) + if n is odd (see [TU] 
Theorem 2.14). In this appendix, we will use the general formula for the second 
derivative in order to reprove this result and to show that the critical points are 
non-degenerate. 

A direct computation of d 2 A n /da 2 (0) using the discriminant works as fol- 
lows: the spectrum of H a is given by A -1 (2 cos a). Near A = A+, we have 
2 + A'(A+)(A n (a) - A+) ~ 2 cos a. This gives A n (a) ~ A+ - a 2 /A'(A+), hence 
the Morse index of A n at a = is if n is odd and 1 is n is even. 
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A direct calculation of the Hessian 



We will denote with a "dot" the derivatives w.r. to a and by a "prime" the 
derivatives w.r. to x. The operator H a is unitarily equivalent to K a = e ~ iax ff e iax 
acting on 1-periodic functions. We have 

K a = H - 2za-f- + a 2 . 

ax 

The derivatives of K a w.r. to a at a = are K = — and K = 2. Applying 
Proposition CD and denoting by (p n a corresponding normalized eigenfunction, we 
get 

A„(0) = 2 + 4i / <p n {x)<f)' n {x)dx . 



Moreover H(j> n (x) = —K(p n = 2i(p' n (x). 

Let us denote by ip the function yi(., 0). Then, using the method of "variation 
of parameters" (i.e. making the Ansatz </> n (x) = Ci(x)i/j(x) + C2(x)(f> n (x) with 
C[ix)ip(x) + C' 2 {x)ct) n {x) = 0), we get 

<fin{x) = -ix(j) n (x) + kip(x) + C(j) n (x) , (5) 

where the constant k is chosen so that (j) n (x) is periodic and C is an arbitrary 
constant which can be fixed by a normalization of 4> n . We can always assume 
that n (O) = n .(l) = by shifting the origin of 1R to some zero of <f> n . Using 
the wronskian, we see that n (l) = n (O). We have to check the derivatives: 
kip'(l) - 1(^(1) + <f>' n (l)) = fc^'(O) - i(f) n (0) or fc^'(l) = z^(0). This gives, using 
Equation ([5]), 

4>n(x) = -ix(f) n (x) + i^j^-ipix) + C(f) n (x) . 

^ (1) 

We get 

A„(0) = 2 + 4z [ [-ix(f> n (x) + ki/){x) + C0 n (x)]</4(x)cfe . 



By integration by parts, we have J* 2x<fi n (x)cj)' n (x)dx = — (f) n (x) 2 dx = — 1. 

Moreover, again by integration by parts, J* ip(x)(f)' n (x)dx = — f i/j'(x)(j) n (x)dx 

and, since the Wronskian tp(f)' n — tp'(f) n is constant and = (j)' n (0), f Q i/j(x)(j)' n (x)dx = 
|^(0). We get 

A n (0) = -20' n (O) 2 M(l) • 

Moreover, it follows from Equation (2.13), page 16 in |10] and the fact that 
4> n = 4>' n (0)y2, that this is exactly — 2/A'(A+). 
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C The case where the eigenfunction vanishes at 
some vertex 



In this Appendix, we take H e Oq and assume that A n = is non-degenerate 
eigenvalue of H with a normalized eigenfunction 0. We have the 

Proposition 2 Let us assume that, for all vertices x satisfying 4>(x) = 0, there 
exists a vertex y ~ x so that (j>{y) ^ 0. Then, for any if; e R x orthogonal to <p, 
there exists a smooth deformation H t £ Oq of H so that <j> = ip. 

It is enough to check that the space of Hip is M, x and to use the first variation 
formulae given in Appendix A. 

Theorem 3 Let us assume that the function <fi vanishes at the unique vertex 
xq. Then, the nullity of the Hessian of the "magnetic variation" of H is at least 
\n + — n_| where n± is the number of vertices x ~ xq so that ±<f)(x) > 0. 

Proof. - 

Choose a smooth variation H t of H so that <f)(xo) = 1. Let v be 
the number of sign changes of <ft away of Xq. Then, for t > small 
enough, the number of sign change of <pt is v + n_ while, for t < 
small enough, it is v + n + . We see from Theorem [1] that the magnetic 
Morse index is v + n_ — in — 1) for t > and v + n + — (n — 1). The 
discontinuity of the Morse index at t = is \n + — n_|. This gives the 
lower bond on the nullity. 

□ 

Corollary 2 // |n+ — n_\ > f3, the eigenvalue is degenerate. 

Let us remark that this lower bound is not always sharp. In the following 
example, we have n + = n_, (3 = 2 and the nullity of the Hessian is 2. 

Example C.l The graph G is made of 2 cycles of length 3 with a common 
vertex. The matrix of H is chosen as follows: 

f 1 1 1 \ 
112 
[H} = - 12 112 
111 
\ 2 1 1 j 

Using the fact that the graph has a symmetry of order 2 exchanging the 2 cycles, 
one can split M. x and the matrix H into the even and odd parts. This allows to 
check that A4 = is non- degenerate. In order to compute the magnetic Hessian, 
we check that it is possible to build a decomposition = d~R x © K which is 

Q orthogonal and with K C kerrf*. It is then easy to check that the magnetic 
Hessian evaluated on K vanishes. 
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D Bipartite graphs 



Let G = (V, E) be a bipartite graph, with V = Y U Z and all edges have a vertex 
in Y, the other in Z. Let U be the involution on MY given by Uf(x) = —f(x) 
if x G Y and Uf(x) = f(x) if x G Z and let 5 be a magnetic field. Then 
UH B U = —H' B with H' G G . So that \\ V \{H B ) = -M h 'b)- An hence ' lt 
follows form the diamagnetic inequality that B — > \\v\(Hb) has a maximum at 
B = 1. And hence the Morse index of the Hessian of — > A|y| (-£/#) at 5 = 1 is 
the dimension of the manifold of magnetic fields namely (3. On the other hand 
the first eigenfunction (f) 1 of H' is everywhere > and the number of sign changes 
of U(pi is \E\. So Berkolaiko's formula for Aiyi gives (|V| — 1) + /3 = \E\. This is 
the Euler formula. 

E Link with the Hessian of the determinant 

Let us assume that we are in the discrete case and the eigenvalue we consider is 
X n = 0. Then we have 

det(H B ) = X n (B)det'(H B ) 

where det^ifg) = F(B) is the product of the eigenvalues Aj for j ^ n. We have 
(-1)«- 1 F(1) > 0. Hence the index of B -> (-l) n ^ 1 det(H B ) is the same as the 
index of B — > \ n (B). 

There is a formula for the characteristic polynomial of a magnetic Laplacian 
on graphs due to Robin Forman [8] and reproved by Kenyon [9] and Burman [3]. 
Using the gauge change / — > f<p n as in my paper gives a Laplace type operator 
whose entries can be of any sign. Forman's formula extends to that case and it 
would be nice to get Berkolaiko's formula form Forman's formula. 
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